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INTRODUCTION 


Optical  communications  through  seawater  can  be  achieved  with  sources 
centered  at  M80  run,  the  blue-green  spectral  region.  Therefore,  crystals 
exhibiting  blue-green  luminescence  are  of  considerable  interest  as  solid 
state  lasers.  Particular  attention  has  been  focused  on  the  Ce+3  d  -*■  f 
luminescent  transition,  which  is  broad  and  appears  to  center  around  520 
nm  in  Y3Ali2°15  (TAG).1'2  The  coordination  in  YAG  has  Ce+3  occupying  the 
center  of  a  cube  of  02-  ions.  In  other  lattices  such  as  SCBO3,  where 
Ce+3  is  coordinated  octahedrally  by  6q2-  ions,  emission  is  typically 
observed  in  the  UV  at  X  <420  nm.2  Blasse  and  Bril2  observed  that  these 
symmetry  sensitive  luminescence  spectra  are  in  keeping  with  sinqsle  con¬ 
siderations  of  crystal  field  theory.  Unfortunately,  Ce+3-YAG  appears  to 
be  an  unsuitable  blue-green  laser  material  because  the  emission  is  at  too 
low  an  energy,  and  because  of  parasitic  excited  state  absorption  processes 
which,  although  not  well-understood,  preclude  laser  action. 

These  accumulated  results  have  led  to  a  search  for  alternative  oxide 
hosts  incorporating  Ce+3  in  a  cubic  8-fold  (or  higher)  coordination.  Two 
new  materials  have  come  to  our  attention.  One  is  0hO2,  in  vhich  Ce+3  would 
have  cubic  site  symmetry,  but  in  which  there  would  be  a  strong  tendency  for 
isovalent  substitution.  However,  our  measurements  on  small  flux -grown 
crystals  of  Ce+3  doped  ThO^  using  pulsed  N2  laser  excitation  (337.1  nm) 
show  a  broad  luminescence  with  Imav  M80  nm.  The  other  host  material, 

LaAlOj ,  shows  an  absorption  onset  of  ^440  nm  when  doped  with  Ce+3.2  This 
would  place  the  d  -*■  f  luminescence  in  approximately  the  correct  spectral 
range.  The  Ce+3  would  have  a  dodecahedral  03- * *  coordination  in  this  crystal. 
We  have  demonstrated  a  flux  synthesis  of  LaA103  (to  stabilize  the  cubic 
structure)  ,  using  La 203,  A12C>3,  and  CeP3  as  reactants  and  2  m/o  BaCC>3  and 
BaF2  as  flux.  The  reaction  is  carried  out  at  11508C  in  air  in  a  Pt  crucible. 
Some  indication  of  blue-green  luminescence  in  Ce+3  doped  La(Al,Sc)03  has 
been  observed  recently.5 

A  critical  factor  in  all  of  these  investigations  is  the  growth  of 
single  crystals  of  these  highly  refractory  host  materials.  Growth  of  Th02 
and  LaA103  can  be  achieved  using  a  suitable  flux,  but  the  crystals  so 
obtained  are  of  irregular  shape  and  must  be  separated  from  the  surrounding 
flux  after  cooling. 


2G.  Blasse  and  A.  Bril,  Appl.  Phys.  Letters  LL,  53  (1967). 

2G.  Blasse  and  A.  Bril,  J.  Chem.  Phys.  47,  5139  (1967). 

3W.  Miniscalco,  J.  Pellegrino  and  W.  Yen,  J.  Appl.  Phys.  49,  6109 

(1978) . 

^W.  Zwicker  et  al.,  ONR  Pinal  Technical  Report,  Contract  N00014-79- 

C-0180 ,  October  1982. 

*V.  Nicolai,  Office  of  Naval  Research,  private  communication. 
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In  order  to  achieve  a  more  controlled  growth  of  these  crystals  by 
the  flux  method,  the  traveling  solvent  zone  (TSZ)  technique  might  be 
employed.  Table  1  lists  seme  oxides  which  have  been  grown  by  this 
method.  A  simple  apparatus  for  crystal  growth  using  this  procedure 
is  shown  in  Figure  1.  The  reactants  are  contained  in  a  boule  with  a 
pointed  tip  to  enhance  crystal  nucleation.  Initially,  the  flux  is 
contained  within  the  tip  of  the  boule.  Above  it  is  a  rod  of  feedstock 
(e.g.,  Ce+3  doped  LaA103)  which  has  been  pressed  and  sintered  into  a 
dense  ceramic  rod.  A  hot  zone  is  defined  by  an  independently  controlled 
heating  element.  As  the  boule  is  lowered  through  the  hot  zone,  dissolu¬ 
tion  of  the  feedstock  occurs  at  the  upper  flux/ceramic  interface  while 
crystallization  of  the  oxide  occurs  at  the  bottom  (existing)  interface. 
The  role  of  "background"  heating  is  uncertain,  but  may  be  necessary  here 
to  maintain  a  hot  zone  temperature  above  ,vl300eC,  if  required;  thus,  it 
is  included  in  Figure  1. 

TABLE  1.  CRYSTALS  OF  OXIDES  GROWN  BY  THE  TRAVELING  SOLVENT  ZONE  METHOD 


Growth 

T  (background)  T (growth  zone)  Rate 
Material  Flux  _ (°C)  _ (°C) _  (mm/hr)  Ref. 


BaTi03 

Ti02 

1290 

1350 

0.1 

6 

CaCC>3 

Li2C03 

ambient 

750 

0.2 

7 

ZnO 

PbF2 

ambient 

900 

0.15 

8 

SrTi03 

PbF2/PbO 

ambient 

1250 

0.05 

9 

(Pb,Sr)Ti03 

PbF2/PbO 

ambient 

1100 

0.05 

9 

PbTi03 

PbF2/PbO 

ambient 

1075 

0.05 

9 

In  this  program,  we  have  evaluated  several  fluxes  for  growth  of 
LaA103,  I*aAll-x,*x03  (M  =  Sc+3,  Ga+3)  and  LaSc03  containing  Ce+3,  and 
have  compared  the  products  to  those  synthesized  in  a  high  temperature 
arc  furnace.  A  TSZ  crystal  growing  furnace  has  been  designed  and  con¬ 
structed,  and  some  preliminary  flux  growth  runs  carried  out.  Modifica¬ 
tions  to  the  crystal  growth  apparatus  have  been  suggested  based  on  the 
results. 


*C.  E.  Turner,  N.  H.  Mason  and  A.  W.  Morris,  J.  Crystal  Growth 
56,  137  (1982). 

7J.  J.  Brissot  and  C.  Belin,  J.  Crystal  Growth  8,  213  (1971) . 
aG.  A.  Wolff  and  H.  E.  LaBelle,  Jr.,  J.  Am.  Ceram.  Soc.  48,  441 
(1965) . 

*B.  DiBenedetto  and  C.  Cronan,  J.  Am.  Ceram.  Soc.  51,  364  (1968). 
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As  a  parallel  project,  a  computer  program  has  been  written  £or  a 
IBM  personal  computer  that  performs  siaple  crystal  field  calculations 
with  the  input  of  various  lattice  parameters.  It  is  anticipated  that 
such  a  program  will  be  useful  in  interpretation  of  host  dependent  Ce+3 
absorption  and  luminescence  spectra.  It  is  reproduced,  with  documenta 
tion,  in  the  Appendix. 
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RESULTS 


Oxide  Syntheses  using  Arc  Furnace  Techniques 

In  order  to  incorporate  Ce+3  into  a  perovskite  lattice,  and  at  the 
same  time  avoid  partial  oxidation  to  Ce+4,  the  synthetic  reactions  must 
be  carried  out  under  neutral  or  reducing  conditions.  One  convenient 
method  for  such  synthesis  is  to  employ  an  arc  furnace.  The  furnace  is 
supplied  commercially  by  Centorr  Assoc.,  Suncook ,  NH,  and  is  essentially 
the  same  design  as  described  by  Reed  and  Pollard10'1 11,  the  originators  of 
the  technique.  The  furnace  body  consists  of  a  thoreated  tungsten  cathode 
centered  above  a  water  cooled  graphite  or  copper  hearth,  which  also  serves 
as  the  anode.  The  arc  is  struck  using  a  conventional  welding  power  supply 
(300  amp,  30V).  The  cathode  and  hearth  are  encased  within  a  quartz  window 
sleeve,  compression  sealed  to  the  top  and  bottom  against  0-ring  gaskets. 

The  cathode  is  attached  to  a  handle  which  protrudes  through  the  top,  held 
by  a  ball  and  socket  joint  which  permits  adjustment  of  the  arc  angle.  Gas 
purge  tubes  are  introduced  into  the  top  and  base,  so  that  the  furnace  can 
be  operated  either  in  vacuo  or  under  a  chosen  atmosphere  -  usually  reducing, 
unless  cathode  and  anode  materials  are  employed  with  high  temperature  oxida¬ 
tion  resistance  (e.g.,  Ir) . 

In  actual  operation,  the  components  of  the  material  to  be  synthesized 
are  introduced  into  a  depression  in  the  hearth.  After  flushing  the  chamber 
and  introducing  the  desired  atmosphere,  the  arc  is  struck  and  directed  into 
the  feedstock.  The  current  is  increased  until  melting  occurs.  Since  the 
hearth  is  water  cooled,  the  molten  charge  "beads  up",  forming  a  thin  skin 
of  solid  material  at  the  interface.  All  of  the  powder  in  the  hearth  may 
be  guided  into  the  bead  by  manipulating  the  arc  position.  The  cold  hearth 
minimizes  contamination  of  the  charge  -  thus,  a  containment  material  such 
as  Cu  may  be  employed,  because  of  its  excellent  thermal  conductivity. 

The  arc  furnace  methods  were  used  for  attempted  syntheses  of  LaA103, 
LaAlo.7Sco.3O3  and  LaSc03,  containing  0.1  to  1.0  mole  %  Ce+3,  added  as 
Ce203.  The  mixtures  were  prepared  stoichiometrically  from  the  component 
oxides.  These  were  ultimately  mixed  and  ground  for  several  hours  using 
an  agate  mortar  in  the  presence  of  MeOH,  to  form  a  jewelers'  paste.  The 
MeOH  was  subsequently  stripped  in  a  vacuum  oven  at  100°C.  The  resulting 
powder  was  cold  pressed  into  pellets  for  introduction  into  the  arc  furnace 
(typically  M.0g  per  run).  All  runs  were  conducted  in  an  Ar  atmosphere. 

All  materials  prepared  in  the  arc  furnace  possessed  a  grey  colora¬ 
tion,  independent  of  the  hearth  used.  If  the  reaction  time  was  too  short. 


1 °T.  B.  Reed,  Mat.  Res.  Bull.  2,  349  (1967). 

11T.  B.  Reed  and  E.  R.  Pollard,  J.  Cryst.  Growth  2,  243  (1968) . 
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there  was  a  tendency  for  the  solidified  bead  to  form  a  "geode,"  with  an 
unreacted  powder  encased  in  a  hard,  amorphous  appearing  shell.  Thus, 
a  reaction  time  of  at  least  10  minutes  was  necessary  to  ensure  complete 
reaction.  'Hie  homogeneity  was  further  promoted  by  stopping  the  fusion, 
turning  over  the  bead,  and  resuming  the  arc  melting  on  what  was  the 
bottom  surface.  (Hiese  techniques  were  practiced  first  by  us  in  the 
laboratory  of  Mr.  Robert  Fahey  ,  at  the  MIT  Lincoln  Laboratories,  before 
being  conducted  on  our  own  apparatus.) 

Samples  prepared  by  the  arc  furnace  technique  were  either  non- 
luminescent  or  only  very  weakly  so,  independent  of  the  Ce  content. 

The  X-ray  diffraction  patterns  of  the  products  are  summarized  in 
Figure  2.  The  LaA103  arc  synthesis  led  unequivocally  to  the  desired 
product,  with  no  X-ray  observable  contaminants  (detection  level  <0.5%) . 
This  was  reproduced  in  several  separate  runs  carried  out  both  at  Lincoln 
Laboratories  and  at  EIC.  Hie  LaSc03*0.1%  Ce  synthesis  led  to  a  product 
with  a  more  complex  X-ray  pattern  than  that  of  the  nearly  cubic  LaAK>3. 

Hie  pattern  contained  the  lines  attributable  to  LaSc03  in  the  orthorhombic 
GdFe03  structure.  Several  impurity  peaks  were  unassigned,  but  appear 
consistent  with  SC2O3.  Definite  assignment  is  difficult  due  to  overlap 
between  most  of  the  major  LaSc03  and  SC2O3  peaks.  La203  was  definitely 
not  present  in  the  product. 

Syntheses  of  LaAlo.7SCQ.3O3  were  carried  out  with  1%  and  0.1%  Ce+3 
added.  Hie  X-ray  patterns  of  the  products  are  of  a  simple  cubic  structure. 
As  seen  in  Table  2,  comparison  with  LaA103  indicates  that  the  mixed  oxide 


TABLE  2.  X-RAY  DIFFRACTION  PARAMETERS  MEASURED  FOR 
ARC  FURNACE  SYNTHESIZED  LaAl03  AND  LaAlo.7SCQ.3O3 


LaA103 

d(A) 

LaAlo.7Sco.3O3 

d  (A) 

hkl 

3.71 

3.86 

012 

2.64 

2.70 

110 

2.17 

2.21 

202,  006 

1.88 

1.91 

024 

1.68 

1.71 

122,  116 

1.54 

1.57 

300,  214, 

1.33 

1.36 

220,  208 

has  a  more  expanded  lattice,  as  would  be  expected  for  the  substitution  of 
Sc+3  (r  =  0.81A)  for  Al+3  (r  *  0.50A).  Again,  the  Ce+3  doped  material 
was  only  very  weakly  fluorescent,  yielding  a  luminescence  spectrum  similar 


12T.  B.  Reed,  R.  E.  Fahey  and  J.  M.  Honig,  Mat.  Res.  Bull.  2,  561 
(1967) . 


7 


X-ray  diffraction  patterns  of  La  oxides  prepared  in  an  arc  furnace  under  Ar 
from  mixtures  of  the  component  oxides.  The  LaAl.03  spectrum  is  in  good  agree 
ment  with  the  literature. 


to  that  reported  earlier6 7 8 9 * * * 13  (Lav  =  450  nm) .  However,  since  the  products 
were  grey,  they  apparently  contained  considerable  concentrations  of  C  or 
reduced  impurity  that  could  lead  to  luminescence  quenching.  From  a 
structural  viewpoint,  the  mixed  La(Sc,Al)03  matrix  would  appear  to  be 
an  interesting  lattice  for  accommodating  Ce+3  in  a  cubic  environment. 

Oxide  Syntheses  from  Fluxes 

For  TSZ  crystal  growth,  it  was  first  necessary  for  us  to  survey 
several  potential  flux  solvents  for  synthesis  of  LaA103  and  La(Al,Sc)03- 
These  experiments  were  conducted  in  a  Lindburg  muffle  furnace  with  a 
controlled  atmosphere  retort.  The  temperature  was  lowered  at  a  pre¬ 
determined  rate  by  a  clock  motor  driven  gear  mechanism  attached  to  the 
rheostat  of  the  temperature  controller.  Syntheses  were  conducted  in 
20  ml  Pt  crucibles  with  mechanically  fitted  tops  (Engelhard  Corp.) . 

In  general,  fluxes  were  selected  which  were  molten  below  'VL400°C, 
the  upper  temperature  limit  of  our  furnace  apparatus.  This  is  also  the 
probable  upper  temperature  limit  that  the  hot  zone  of  a  boule  lowered 
through  our  TSZ  furnace  could  be  expected  to  attain  homogeneously  (see 
below) .  It  is  also  of  interest  to  seek  low  temperature  growth  of  LaAlC>3 
single  crystals,  in  order  to  avoid  the  twin-producing  cubic  to  tetragonal 
phase  transition  at  450°C.11* 

Representative  fluxes  and  conditions  of  syntheses  are  summarized  in 
Table  3.  In  all  cases,  the  starting  materials  consisted  of  mixtures  of 
the  component  oxides  (e.g.,  1^303,  AI2O3  and  06303)  plus  a  molar  excess 
of  the  flux.  The  ratio  of  feedstock  material  to  flux  was  chosen  accord¬ 
ing  to  various  literature  sources  for  flux  growth  of  perovskite  crys¬ 
tals.6  9 

The  BaF2/BaC03  flux  (not  shown  in  Table  3)  did  not  melt  at  the 
reaction  temperature.  This  material  does,  however,  promote  the  solid 
state  reaction  of  1*6303  and  AI2O3  to  yield  LaA103.  We  have  previously 
reported  the  synthesis  of  blue-green  luminescent  Ce+3-doped  LaA103 
using  this  technique.13 

The  two  oxide  fluxes,  ^38407  and  481303*7305,  did  form  a  molten 
flux.  The  components  of  the  desired  compound  (e.g.,  83303,  AI2O3,  etc.) 
were  mixed  with  the  flux  by  intimate  grinding  with  an  agate  mo-tar  in  a 
MeOH  paste.  The  mole  percent  of  flux  in  these  syntheses  was  between 


6Tumer  et  al. ,  loc.  cit. 

7Brissot  and  Belin,  loc.  cit. 

8Wolff  and  LaBelle,  loc.  cit. 

9DiBenedetto  and  Cronan,  loc.  cit. 

1SR.  D.  Rauh  and  A.  N.  Scoville,  ONR  Final  Technical  Report,  Contract 

N00014-82-C-0142,  July  1982. 

16J.  P.  Coutures,  J.  M.  Badi,  R.  Berjoan,  J.  Coutures,  R.  Flamand  and 

A.  Rouanet,  High  Temp.  Science  13,  331  (1980). 
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40  and  60%.  The  oxide-flux  mixtures  were  allowed  to  equilibrate  at  the 
upper  temperature  (1200  to  1400°C)  for  several  hours  before  cooling  at 
3-20°C/hr  to  a  temperature  near  the  flux  solid  function  (750-900°C) . 

Then,  the  fluxes  were  either  rapidly  quenched  or  poured  off  any  product 
as  a  melt.  The  reactions  with  these  oxide  fluxes  were  all  conducted  in 
an  air  atmosphere. 

X-ray  diffraction  patterns  of  the  products  of  the  oxide  flux  reac¬ 
tions  indicated  that  AI2O3  was  a  major  product  component,  and  that  LaAlC>3 
and  LaSc03  were  minor  components,  if  present  at  all.  For  the  most  part, 
the  products  of  the  flux  reactions  had  different  X-ray  diffraction  patterns 
from  the  arc-furnace  prepared  materials.  Although  there  were  some  similar¬ 
ities  between  the  X-ray  patterns  of  the  4Bi203*V205  flux  products  and  that 
of  LaA103,  several  peaks  are  missing  or  shifted,  so  that  positive  identifica¬ 
tion  is  questionable.  The  X-ray  pattern  of  this  product  is  shown  in  Figure  3. 

KF  also  remained  molten  throughout  the  reaction  temperature  range. 

The  X-ray  pattern  of  the  product  of  the  reaction  of  La2°3  and  AI2O3  in  KF 
was  identical  to  that  produced  in  the  arc  furnace  synthesis,  and  is  iden¬ 
tified  as  LaA103.  The  product  was  a  white  powder  and  did  not  fluoresce 
despite  the  fact  that  the  reaction  was  carried  out  under  Ar.  The  sample 
of  LaAlo.7Sco.3O3,  1%  Ce  prepared  in  a  KF  flux  had  a  complicated  X-ray 
spectrum,  but  only  one  dominant  line  at  d  *  2.84 &.  The  spectrum,  shown 
in  Figure  4,  is  similar  to  that  of  the  garnet  La3Sc2Ga3012 >  which  can  be 
prepared  from  the  solid  state  reaction  of  La203»  SC2O3  and  Ga203.  This 
type  of  phase  is  not  unlikely  in  the  present  case,  since  Al+3  and  Ga+3 
often  substitute  isostructurally . 

Several  experiments  were  conducted  on  growth  of  LaAli_xGai_x03 
using  the  fluoride  and  oxide  fluxes.  However,  weighing  the  crucible 
before  and  after  the  syntheses  indicated  that  the  Ga203  had  been  lost 
to  volatilization.  Hence,  no  further  experiments  were  performed  on  this 
system. 

Traveling  Solvent  Zone  Experiments 

1.  Furnace  Design  and  Assembly .  The  traveling  solvent  zone 
furnace  takes  advantage  of  the  temperature  dependent  solubility  of  re¬ 
fractory  materials  in  various  fluxes.  A  cylindrical  boule  with  a  coni¬ 
cally  tapered  base  is  filled  from  bottom  to  top  with  a  seed  crystal, 
a  layer  of  flux,  and  a  dense  column  of  feed  (oxide)  material.  The  whole 
boule  is  heated  up  to  just  below  the  melting  point  of  the  flux,  then  a 
hot  ring  placed  exactly  at  the  level  of  the  flux  layer  heats  this  thin 
zone  of  the  boule  to  a  higher  temperature  so  as  to  melt  the  flux.  The 
boule  is  then  dropped  slowly  (e.g.,  <1  cm/day)  down  through  the  hot  ring. 

The  molten  layer  of  flux  dissolves  the  feed  material  above  and  deposits 
a  single  crystal  below. 

The  basic  components  of  the  furnace  are  the  frame,  the  heating 
column,  the  clock  drive,  and  the  power  supplies.  The  frame  is  a  simple 
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13”  x  24"  x  24"  angle-iron  assembly  with  plenty  of  holes  and  crossbars 
for  attaching  parts.  The  other  components  are  described  in  detail. 

The  heating  column  is  a  three-zone  tube  furnace.  The  identical 
outer  zones  consist  of  Kanthal  REH  B  7-30  tube  elements  wrapped  in  a 
2"  blanket  of  Kaowool  ceramic  insulation  and  secured  by  ring  clamps 
(see  photograph,  Figure  5) .  The  middle  zone  is  a  1"  diameter  Kanthal 
K  3/6  "Superelement,"  a  thin  MoSi2  ring  capable  of  attaining  temperatures 
in  excess  of  1500°C.  This  element  was  fitted  with  a  ceramic  retaining 
block  which  served  multiple  functions:  1)  to  retain  and  position  the 
fragile  ring  element;  2)  to  expose  only  a  narrow  slit  of  the  ring  element 
to  the  furnace  center?  3)  to  position  the  outer  zones  of  the  furnace  con¬ 
centrically  with  each  other  and  the  ring.  This  last  function  was  achieved 
by  machining  helical  grooves  into  the  holder  which  match  the  ends  of  the 
Kanthal  tube  elements.  The  retaining  block  material  was  a  soft,  low 
density  refractory,  Salir,  manufactured  by  Zircar  Corp.  The  ring  element 
and  retaining  block  are  shown  in  the  accompanying  photographs  (Figure  6) . 

Alumina  tubes  were  passed  through  the  Kanthal  elements  and  positioned 
by  a  groove  in  the  ring  element  retainer.  These  protect  the  boule  from 
electrical  hazards  and  thermal  discontinuities  presented  by  the  Kanthal 
heating  coils.  A  stainless  steel  adaptor  was  attached  to  the  bottom  of 
the  column  to  allow  passage  of  argon  through  the  furnace. 

Figure  7  is  a  drawing  of  the  hot  ring  and  retainer,  the  upper  alumina 
tube,  and  the  boule.  The  boule,  shown  in  its  loading  position  at  the  top 
of  the  furnace,  is  a  vitreous  carbon  cylinder  tapered  to  a  point  at  the 
bottom  to  allow  insertion  of  a  seed  crystal.  A  carriage  made  of  stainless 
steel  washers  and  oxidation  resistant  metal  wires  (e.g. ,  W,  Ni,  or  Mo) 
attach  the  boule  to  the  clock-drive.  The  clock  is  a  1  or  1/12  rph  motor 
with  a  gear  train  which  allows  selection  of  a  variety  of  dropping  rates 
in  the  range  of  1  to  50  mm  per  day. 

Figure  8  is  a  schematic  of  the  power  supply  to  the  outer  and  center 
zones  of  the  furnace.  The  Kanthal  tube  elements  operate  at  a  nominal 
24V/58A,  and  were  driven  by  high-current  transformers  powered  by  small, 

10A  variacs  (20A  variacs  will  be  needed  if  the  elements  are  to  be  operated 
near  full  power) .  The  ring  element  requires  a  nominal  6V/40A  and  was 
powered  by  an  externally  programmable  Sorenson  SRL  40-50  DC  power  supply 
regulated  by  an  LFE  temperature  controller.  These  instruments  were  com¬ 
patible  with  no  fuss,  and  excellent  temperature  regulation  was  achieved. 

Temperature  profiles  of  the  furnace  were  obtained  by  dropping  a 
thermocouple  down  the  center  tube  and  measuring  temperature  as  a  function 
of  distance  from  the  top  of  the  furnace.  The  variables  which  determine 
the  shape  of  this  profile  are  the  power  applied  to  the  outer  zones,  the 
temperature  setting  of  the  hot  ring,  and  the  rate  of  flow  of  argon  or 
other  gas  through  the  furnace.  It  is  expected  that  the  presence  of  the 
boule  in  the  furnace  will  radically  alter  the  profile  shape.  It  is 
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Figure  5.  Photograph  showing  disassembled  TSZ  crystal 

growing  furnace,  including  insulation  wrapped 
2-zone  furnace  and  supports,  the  central 
ceramic  guide  tubes  with  inert  atmosphere 
purge  adaptor,  and  the  central  MoSi2  ring 
heating  element  with  ceramic  retainer. 
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desirable  to  have  a  sharply  defined  temperature  peak,  with  the  high  tempera¬ 
ture  zone  (i.e.,  the  region  over  which  T  is  sufficient  to  melt  the  flux) 
occurring  over  a  vertical  distance  of  <1  cm.  Figure  9  depicts  the  best 
profile  obtained  to  date  from  our  furnace.  Supposing  we  were  using  a 
potassium  fluoride  flux  <mp.  *  850°C) ,  the  high  temperature  zone  occurs 
over  a  vertical  distance  of  1  inch.  This  profile  was  obtained  with  the 
outer-zone  variacs  set  to  15V  (and  thoroughly  equilibrated) ,  the  ring 
temperature  at  1250°C,  and  no  Ar  flow. 

2.  TS2  Flux  Growth  Experiments.  Several  trial  runs  have  been  made 
in  our  furance.  None  of  these  runs  were  successful  in  producing  single 
crystals. 

The  first  series  of  runs  used  an  open  boule  containing  an  LaA103 
oxide  mix  with  10  m/o  KF  flux.  Very  soon  after  turning  on  the  hot  ring, 
the  Mo  filament  holding  the  boule  gave,  apparently  having  been  corroded 
by  flux  vapor. 

In  the  second  trial,  the  KF  flux  was  not  mixed  with  the  oxide,  but 
was  melted  into  a  plug  at  the  bottom  of  the  boule.  A  platinum  foil  cap 
was  wired  to  the  top  of  the  boule,  and  a  tungsten  filament  replaced  the 
molybdenum  of  run  1.  Approximately  4  hours  after  turning  on  the  hot  ring, 
the  boule  fell  to  the  bottom  of  the  furnace.  This  time  the  boule  itself 
had  oxidized,  indicating  that  measures  must  be  taken  to  seal  the  furnace 
interior  from  the  atmosphere.  The  material  at  the  tip  of  the  boule  con¬ 
tained  some  small,  nonluminescent  crystals,  <1  mm  to  an  edge.  These  were 
analyzed  as  LaAl03  by  x-ray  diffraction,  it  soon  became  obvious  that 
volatility  of  the  fluoride  flux  was  a  major  problem.  Time  permitted  us 
only  one  attempted  solution  at  operating  the  growth  experiment  in  a  closed 
system.  In  this  experiment,  a  cylindrical  nickel  crucible  containing  a 
plug  of  KF  flux,  a  stoichiometric  (LaA103)  oxide  mix,  and  capped  with  a 
nickel  cap,  was  vacuum-sealed  inside  a  quartz  ampoule.  This  ampot ie  was 
then  dropped  down  the  furnace.  Approximately  12  hours  after  turning  on 
the  hot  ring,  the  ampoule  exploded,  apparently  due  to  a  buildup  of  internal 
pressure  and  a  weakening  of  the  quartz  walls  by  fluoride  vapor. 

The  furnace  has  been  repaired  and  supplied  with  a  fireproof  ceramic 
base,  and  is  ready  for  further  use. 
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Figure  9.  Typical  temperature  profile  for  crystal  growth  furnace. 
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CONCLUSIONS 


Although  the  crystal  growth  experiments  did  not  yield  crystals  suit 
able  for  optical  evaluation,  the  compound  synthesis,  furnace  construction 
and  trial  TSZ  runs  did  provide  a  basis  for  further  research.  We  draw  the 
following  conclusions  from  these  experiments: 

•  A  single  phase  cubic  mixed  oxide  of  formula  LaAlxSci_x03 
can  be  produced  by  arc  melting  of  stoichiometric  mixtures 
of  the  component  oxides.  The  lattice  dimensions  are 
expanded  compared  to  LaA103 . 

•  Although  synthesis  of  LaA103  and  LaAl]_-xScx03  from  oxide 
fluxes  was  unsuccessful  here,  it  has  been  reported  else¬ 
where.15  However,  oxidizing  fluxes  are  not  likely  to  be 
compatible  with  Ce203»  which  is  easily  oxidized  itself. 

Hence,  we  chose  to  emphasize  the  neutral  flux,  KF,  which 
has  been  used  to  grow  other  perovskites. 

•  LaAlC>3  was  produced  successfully  in  KF  between  1250  and 
750°C.  The  mixed  LaAli-xScx03  may  have  produced  a  garnet 
phase,  which  in  itself  would  be  a  candidate  cubic  crystal 
host  for  Ce+3. 

•  The  only  luminescent  Ce-doped  oxides  produced  so  far  in 
this  laboratory  have  employed  CeF3  as  a  dopant  source. 

It  may  be  that  Ce2C>3,  which  cannot  be  obtained  in  high 
purity,  is  so  sensitive  to  oxidation  that,  even  under 
neutral  conditions,  residual  oxygen  in  the  growth  appa¬ 
ratus  is  sufficient  to  convert  it  to  CeC>2. 

•  The  failure  of  the  TSZ  growth  was  due  to  KF  flux  vola¬ 
tility  and  reactivity,  as  well  as  apparent  air  leaks 
into  the  growth  chamber.  We  would  make  the  following 
modifications  to  the  growth  apparatus  in  order  to  prevent 
malfunction: 

1)  The  charge  should  be  hermetically  sealed  in  a  KF- 
resistant  thin  walled  boule,  e.g.,  Pt  or  perhaps 
Ni  or  Zr.  This  could  be  done  in  an  arc  welder 
under  inert  atmosphere. 


15T.  Takeda,  T.  Miyata,  F.  Muramatsu  and  T.  Tomiki,  J.  Electrochem. 
Soc.  127,  438  (1980). 


2)  The  boule  should  be  sealed  into  an  evacuated  quartz 
tube  to  eliminate  the  necessity  of  making  the  growth 
furnace  absolutely  airtight. 

3)  The  use  of  rf  heating  of  the  hot  zone  should  be  con¬ 
sistent,  since  the  thermally  heated  zone  had  a  rather 
broad  temperature  profile. 
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APPENDIX 


During  this  program  we  have  begun  an  examination  of  the  theoretical 
basis  for  interpretation  of  Ce+3  luminescence  in  various  hosts.  We  note 
that  green  luminescence  has  been  well  documented  for  Ce  in  yttrium  alumina 
garnet,  where  0=  coordination  is  eightfold  cubic.  In  addition,  blue-green 
luminescence  of  Ce+3  appears  to  have  been  observed  in  the  cubic  perovskite 
hosts  LaAl03  and  LaAlxSci_x03 ,  with  a  dodecahedral  0=  coordination  in  the 
cubic  “A"  site  of  ABO3  lattice.  However,  the  emission  spectrum  of  Ce+3 
in  the  perovskite  YAIO3  is  in  the  ultraviolet  spectral  region  (Xmax =*360  nm) . 
YAIO3  has  a  distorted  perovskite  structure,  or  D2h  symmetry,  similar  to 
GdFe03  (Figure  4) ,  while  LaA103  is  nearly  cubic.  Thus,  the  Ce+3  fluores¬ 
cence  energy  appears  to  be  very  sensitive  to  rather  small  changes  in  the 
crystalline  environment. 

It  would  be  valuable  to  have  a  ready  means  of  predicting  the  effects 
of  hot  lattice  parameters  and  geometry  on  Ce+3  energy  levels,  as  an  aid 
to  interpretation  of  luminescence  and  absorption  spectra.  To  this  end, 
we  have  written  a  simple  crystal  field  program  for  use  with  the  IBM  personal 
computer  with  128K  of  memory.  The  program  is  written  in  the  IBM  Fortran 
language,  and  is  thus  easily  modified  to  run  on  larger  IBM  computers.  It 
is  reproduced  at  the  end  of  this  appendix,  and  is  fully  documented  so  as 
to  be  self-explanatory.  A  copy  of  the  disk  is  available  from  EIC  Labora¬ 
tories. 

In  operating  the  program  on  a  standard  IBM  PC,  the  calculation  can 
be  quite  slow.  A  small  change  in  the  calculation  of  the  matrix  elements 
translates  to  a  larger  change  in  the  resultant  eigenvalues.  The  conse¬ 
quence  of  this  is  that  in  order  to  calculate  the  energy  levels  to  within 
^10%,  the  matrix  elements  must  be  computed  to  <1*  accuracy.  A  Monte  Carlo 
technique  is  used  to  compute  the  numerical  values  of  the  integrals.  Our 
benchmark  tests  of  the  program  indicate  that  in  order  to  achieve  such 
precision,  a  total  of  1.5  x  106  Monte  Carlo  steps  are  required  for  tenfold 
coordination.  Other  integration  algorithms  may  be  possible,  but  this  is 
obviously  the  "rate  determining  step."  A  recent  development  which  can  be 
used  with  the  IBM  PC  is  the  INTEL  8087  coprocessor  chip  for  high  speed 
floating  point  numerical  operations.  This  device  should  be  capable  of 
enhancing  the  speed  of  the  integration  by  more  than  a  factor  of  10. 


24 


c 

c 

c 


cat  eic.f 


MAIN  PROGRAM 


C  CRYSTAL  FIELD  CALCULATION  FOR  CE  +  3  ION 

C 

C - 


c 

c 

C  PURPOSE? THIS  PROGRAM  PERFORMS  DEGENERATE  PERTURBATION  THEORY 
C  CALCULATIONS  FOR  THE  5-D  LEVELS  OF  ANY  ION  IN  THE 

C  ELECTROSTATIC  FIELD  OF  OF  ANY  OTHER  IONS.  THE  POTENTIAL 

C  OF  INTERACTION  MUST  BE  SUPPLIED  IN  AN  FUNCTION  POT(X,Y.Z) 

C  THIS  POTENTIAL  FUNCTION  CAN  BE  COMPILED  SEPARATELY  AND 

C  LINKED  TO  THE  OBJECT  CODE  OF  THIS  PROGRAM  USINb  THE  IBM 

C  LINKER.  THE  POTENTIAL  FUNCTION  MUST  BE  EXPANDED  IN  A 

C  TAYLOR-SERIES  AS  DESCRIBED  IN  TINKHAMS  BOOK  'GROUP 

C  THEORY  AND  QUANTUM  MECHANICS  *  PS  68.  ONLY  THE  FIRST 

C  TERM  WITH  NONSPHERICAL  SYMMETRY  NEED  BE  KEPT. 

C 

C - 


C 

C  METHOD?  THE  PROGRAMS  OPERATES  IN  THE  FOLLOWING  SEQUENCE.  FIRST 
C  THE  COMPUTATION  PARAMETERS  ARE  ESTABLISHED.  THE  NEXT 

C  STEP  IS  TO  COMPUTE  THE  MATRIX  ELEMTENTS  <5DIVI5D>. 

C  SINCE  THE  MATRIX  IS  BOTH  REAL  AND  SYMMETRIC  >  ONLY  15 

C  MATRIX  ELEMENTS  ARE  COMPUTED.  THE  3-DIMENSIONAL 

C  INTEGRALS  ARE  COMPUTED  USING  A  MONTE-CARLO  INTEGRATION 

C  METHOD.  TO  ACHEIVE  AN  ACCURACY  OF  57.  APPROXIMATELY 

C  5000  POINTS  ARE  REQUIRED.  THE  RANDOM  NUMBERS  ARE  READ 

C  INTO  THE  ARRAYS  X,Y»Z.  THIS  WAS  DONE  TO  IMPROVE  THE 

C  SPEED  OF  THE  CALCULAT I ON . THE  FINAL  STEP  IS  TO  COMPUTE 

C  THE  EIGENVALUES.  THE  SUBROUTINE  EIGEN  ACCOMPLISHES  THIS 

C  USING  THE  DIAGONALIZATION  METHOD  ORIGINATED  BY 

C  JACOBI  AND  ADAPTED  BY  VON  NEUMANN  FOR  LARGE  COMPUTERS. 

C 

C - 


C 

C  USAGE:  THE  FOLLOWING  SUBROUTINES  ARE  USED 

C  EIGEN:  COMPUTES  EIGENVALUES  OF  A  REAL  SYMMETRIC  MATRIX 

C  INTEG:  MONTE-CARLO  INTEGRATION  ROUTINE 

c  chang:cqnverts  matrix  m  to  storage  mode=i  matrix  a 

c 

C  THE  FOLLOWING  FUNCTION  ARE  USED 

C  F5DXY»F5DYZ»F5DZXfF5DZ2»F5D22 

C  G5DXY»G5DYZfG5DZX»G5DZ2»G5D22 

C 


C 

C - 

C 

C  REQUIREMENTS 

C 

C 

r - 


AN  ADDITIONAL  FUNCTION  POT(X»Y»Z)  MUST  BE  SUPPLIED 
SEE  DOCUMENTATION  FOR  DETAILS. 
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MAIN  PROGRAM 


STORAGE  J  2 

PROGRAM  CRYST 

EXTERNAL  F5DXY . F5DYZ . F5DZX . F5DZ2 » F5D22 

__  -  -EXTERNAL -GSUXYLi-GSIltZuGSItZX.  G5EZ2  *-G5Il22 _ 

REAL*4  M<5»5> »A(15) »R(15) 

COMMON  /RNDDAT /  X ( 5000 )» Y ( 5000 >* Z ( 5000  ) 

COMMON  /PARAM/ZE.ZQtBD 


INPUT  CALCULATION  PARAMETERS 


5  WRITE(#*10) 

10  FORMAT ( '  ENTER  0  TO  RUN  WITH  INTERNAL 
INTERNAL  PARAMETERS  ARE*. 

CHARGE  ON  CE+=2.0 


PARAMETERS' ) 


20 


30 


CHARGE  ON  0  =-2.0 
LATTICE  PARAMETER=2 . 0 


EtOHRS 


40 

50 

60 

70 

80 

90 

100 

99 


WRITE(#*20) 

FORMAT ( '  ENTER  1  TO  SPECIFY  NEW  PARAMETERS') 
READ(#*30) ICHOIC 
FORMATt II ) 

IF  ( ICHOIC . EQ < 1 ) GOTO  40 

IF( ICHOIC. NE.O)GOTO  5 

BD=2 . 0 

ZE  =  2 . 0 

Z0=2 . 0 

GOTO  99 

WRITE(#*50) 

FORMAT ( '  INPUT  EFFECTIVE  NUCLEAR  CHARGE  F4.2') 
READ ( * *  60 ) ZE 
FORMAT (F4.2) 

WR I TE  ( #  »  70 ) 

FORMAT < '  INPUT  CHARGE  ON  OXYGEN  F4.2') 

READ(*»80)Z0 

FORMAT ( F4 . 2 ) 

WRITE (#*90) 

FORMAT ( '  LATTICE  PARAMETER  F4.2') 

READC#* 100)BD 
FORMAT ( F4 . 2 ) 

CONTINUE 


COMPUTE  MATRIX  ELEMENTS 


WRITE ( #  *  23 ) 

23  FORMAT ( '  HOW  MANY  MONTE-CARLO  POINTS  (14)  LESS  THAN  5000' 
READ(#*24)NSTEP 

24  FORMAT ( 14 ) 

WRITE(#» 333) 

333  FORMAT ( '  READING  IN  MONTE-CARLO  DATA') 
QP£N(1*FILE='RANDQM' ) 
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C  READ  IN  THE  RANDOM  NUMBERS  AND  SCALE  THEM  TO  SPAN  -20  TO  +20 
DO  444  1  =  1  f  NSTEP 
READ < 1.555) X( I ) »Y(I ) »Z(I) 

X(I)«X(I)*40. 0-20.0 
Y(I)=Y(I)*40. 0-20.0 
Z(I)=Z(I>*40. 0-20.0 
555  FORMAT ( 3F8 . 7 ) 

444  CONTINUE 

CLOSE (1fSTATUS=' KEEP  '  ) 

0PEN(2fFILE='RESULTS'  ) 


COMPUTE  MATRIX  ELEMENTS 

EACH  CALL  OF  INTEG  COMPUTES  ONE  MATRIX  ELEMENT 
WRITE  <  *  f 158 ) 

CALL  INTEG<F5DXYfG5DXYfM< 1 f 1 )  f NSTEP) 
WRITE<*» 158) 

CALL  INTEG(F5DXYfF5DYZfM(1f2>  fNSTEP) 

WRITE ( *  > 158 ) 

CALL  INTEG(F5DYZi G5DYZfM<2f 2> , NSTEP) 
++RITEO#  158)  - 

CALL  INTEG<F5DXYfF5DZXfM(1f3) fNSTEP) 

WRITE ( *  r 158 ) 

CALL  INTEG<F5DYZfF5DZXfM(2f3) fNSTEP) 

WRITE ( *  f 158 ) 

CALL  INTEG ( F5DZX fG5DZXfM(3f3)f NSTEP) 

WRITE 158 ) 

CALL  INTEG ( F5DXY , F5DZ2 » M < 1.4) t NSTEP) 

WRITE ( #  f 158  > 

CALL  INTEG(F5DYZ»F5DZ2»M(2»4) » NSTEP ) 

WRITE ( * » 158 ) 

CALL  INTEG ( F5DZX f  F5DZ2 >M<3»4) » NSTEP) 
WRITE(*,158) 

CALL  INTEG  <  F5DZ2  f G5DZ2 fM(4f4)f NSTEP) 

WRITE  <  *  *  158 ) 

CALL  INTEG(F5DXYfF5D22fM<1f5) fNSTEP) 

WRITE  <  * » 158 ) 

CALL  INTEG(F5DYZfF5D22fM(2f5)fNSTEP> 
WRITE(*f 158) 

CALL  INTEG(F5DZX»F5D22»M(3?5) fNSTEP ) 
WRITE<*» 158) 

CALL  INTEG<F5DZ2fF5D22fM< 4  f5) fNSTEP) 
WRITE(*» 158) 

CALL  INTEG ( F5D22 » G5D22 f  M  ( 5  f  5  )  fNSTEP) 

158  FORMAT  (  '  COMPUTING  MATRIX  ELEMENT') 

SYMMETIZE  MATRIX  M 


150 


DO  150  J= 1 r  5 
DO  150  I  =  1fJ 
M  <  J » I )  *M  ( I » J ) 

COMPUTE  EIGENVALUES 
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J 


WRITE!*) 159) 

159  FORMAT!'  COMPUTING  EIGENVALUES') 

N  =  5 
MV=  1 

C  FIRST  CONVERTS  MATRIX  M  TO  STORAGE  MODE  1  MATRIX  A 
CALL  CHANG  !  M  >  A  ) 

CALL  EIGEN ( A  >  R  »  N  » MV ) 

C 

WRITE ! *  >  63 ) 

WRITE ! 2  >  63 ) 

63  FORMAT ( '  THE  CALCULATION  PARAMETERS  ARE . ') 

WRITE ! *  >  64 ) ZO 

WRITE! 2 » 64 ) ZO 

64  FORMAT!/'  CHARGE  ON  OXYGEN  '»F4.2) 

UR  ITE !  *  >  65  >  ZE 

WRITE ! 2 » 65 ) ZE 

65  FORMAT!'  EFFECTIVE  NUCLEAR  CHARGE  '»F4.2) 

UR ITE ! *  >  66 ) BD 

WRITE ! 2  >  66 ) BD 

66  FORMAT!'  LATTICE  PARAMETER  '>F4.2) 

WRITE ! *  >  20 1 ) 

WRITE ! 2 »  201 ) 

201  FORMAT!'  THE  EIGENVALUES  ARE . ') 

URITE!*>645)A!1) >  A  !  3  > »A!6)»A!10)»A!15) 

WRITE! 2 >645) A! 1) >A!3)>A!6)>A!10)>A!15) 

645  FORMAT ! 5E14 . 6 ) 

CLQSE!2>STATUS='KEEP' ) 

END 

C 

C  INTEGRATION  SUBROUTINE 
C 

C  THIS  SUBROUTINE  COMPUTES  THE  3-D  INTEGRAL  OF  F1*P0T*F2 
C  WHERE  FI  AND  FZ  ARE  EXTERNAL  FUNCTIONS  REPRESENTING  THE 
C  ANY  TWO  OF  THE  5D  ORBITALS  AND  POT  IS  THE  POTENTIAL 
C  DUE  TO  THE  COORDINATED  OXYGEN  IONS. 

C  THE  INTEGRATION  METHOD  IS  MONTE-CARLO.  THE  RANDOM  NUMBERS 
C  ARE  STORED  IN  THE  ARRAYS  X>Y>Z 
C 

SUBROUTINE  INTEG ! F 1 > F2 > VALU > NSTEP ) 

EXTERNAL  F1>F2 

COMMON  /RNDDAT /  X ! 5000 )> Y ! 5000 )> Z ( 5000 ) 

COMMON  /PARAM/ZE  >  ZO  >  BD 

CQNST=SQRT !3.00)/1536.0/SQRT !3. 14159265) 

C0NST=C0NST*ZE**3.5 

VALU=0 • 0 

RANGE=40 .0 

DO  20  1  =  1 >  NSTEP 

VALU=VALU+F1!X! I) >Y! I ) » Z ! I ) > ZE ) *F2 ! X! I ) > Y ! I ) , Z ! I )  > ZE > 
9*P0T!X!I)»Y!I)>Z!I)) 

20  CONTINUE 

VALU=VALU*C0NST*C0NST*40.00/FL0AT! NSTEP) 
URITE!*>245)VALU 
URITE!2>245)VALU 
245  FORMAT ! E 14 , 6 ) 

END 
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SUBROUTINE  CHANG (  M »  A  ) 

C  CONVERT  THE  MATRIX  M  TO  MATRIX  A  (STORAGE  M0DE=1>  FOR  USE 
C  BY  EIGEN. 

REAL44  M<5» 5) » A(15) 

K*1 

00  10  J*1 1 5 
00  10  I»1>J 
A(K)*M(I> J) 

10  K*K+1 

END 

5D  ORBITAL  FUNCTIONS 
SDXY  ORBITAL 

FUNCTION  F5DXY. (XrYtZrZE) 

R2«X*X+Y*Y+Z*Z 
R*SQRT ( R2 ) 

F5DXY»<6.0-ZE*R/2.0>*2.0*X*Y*EXP<-1.0*ZE*R/4.0> 

END 

C  SDYZ  ORBITAL 

FUNCTION  FSDYZ (XiYiZiZE) 

R2-X*X+Y*Y+Z*Z 

R«SQRT(R2) 

F5DYZ«(6.0-ZE*R/2.0)*2.0*Y*Z*EXP<-1.0*ZE*R/4.0> 

END 

C  SDZX  ORBITAL 

FUNCTION  F5DZX ( X » Y  » Z » ZE ) 

R2»X*X+Y*Y+Z*Z 
R*SQRT (R2) 

F5DZX«<6.0-ZE*R/2.0)t2.0*Z*X*EXP<-1.0*ZE*R/4.0> 

END 

C  5DZXZ  ORBITAL 

FUNCTION  F5DZ2(X*Y  »Z»ZE) 

R2»X*X+Y*Y+Z*Z 

R-SQRT(R2> 

T-(3.0*ZtZ-R2>*C6.0-ZE*R/2.0>/SQRTC3.0> 

F5DZ2»T*EXP ( -1 . 0*ZE*R/4 . 0 ) 

END 

C  3DX*X-Y*Y  ORBITAL 

FUNCTION  F3D22 ( X »  Y  t  Z  r  ZE ) 

R2»X*X+Y*Y+Z*Z 
R=»SQRT  (R2 ) 

- - T»<X*X“Y*Y)*(6.0-ZE*R/2.0<M-  -  - 

F5D22»T*EXP ( -1 . 0*ZE*R/4 . 0 ) 

END 

C  SDXY  ORBITAL 

FUNCTION  GSDXY ( X  »  Y  »  Z »  ZE ) 

R2»X*X+Y*Y+Z*Z 

R»SQRT(R2) 

G5BXY»(6.0-ZE*R/2.0>*2.0*X*Y*EXP(-1.0*ZE*R/4.0> 

.END 
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C  5DYZ  ORBITAL 

FUNCTION  G5DYZ (  X  »  Y  »  Z  »  ZE ) 

R2*X*X+Y*Y+Z*Z 

R=SQRT(R2) 

G5DYZ«(6.0-ZE*R/2.0)*2.0*Y*Z*EXP(-1 .0*ZE*R/4.0> 
END 

C  SDZX  ORBITAL 

FUNTION  B5DZX(X»Y» Z.ZE) 

R2*X*X+Y*Y+Z*Z 

R=SQRT(R2) 

C  R*SQRT (  R2 ) 

T*(3.0*Z*Z-R2)*( A.0-ZE*R/2.0)/SQRT(3.0) 
G5DZ2»T*EXP(-1.0*ZE*R/4.0> 

END 

C  5DX*X-Y*Y  ORBITAL 

FUNCTION  G5D22 ( X  » Y  » Z » ZE ) 

R2«X*X+Y*Y+Z*Z 
R=SQRT ( R2 ) 

T*  (  X*X-Y*Y ) *  <  6  *  0-ZE#R/2 . 00 ) 

G5D22*T*EXP  ( -1 .0*ZE*R/4 . 0 ) 

END 


SUBROUTINE  EIGEN 
PURPOSE 

COMPUTE  EIGENVALUES  AND  EIGENVECTORS  OF  A  REAL  SYMMETRIC 
MATRIX 

USAGE 

CALL  EIGEN(AfR'NfMU) 

DESCRIPTION  OF  PARAMETERS 

A  -  ORIGINAL  MATRIX  (SYMMETRIC) »  DESTROYED  IN  COMPUTATION. 
RESULTANT  EIGENVALUES  ARE  DEVELOPED  IN  DIAGONAL  OF 
MATRIX  A  IN  DESCENDING  ORDER. 

R  -  RESULTANT  MATRIX  OF  EIGENVECTORS  (STORED  COLUMNWISE  t 
IN  SAME  SEQUENCE  AS  EIGENVALUES) 

N  -  ORDER  OF  MATRICES  A  AND  R 
MV-  INPUT  CODE 

0  COMPUTE  EIGENVALUES  AND  EIGENVECTORS 
1  COMPUTE  EIGENVALUES  ONLY  (R  NEED  NOT  BE 

DIMENSIONED  BUT  MUST  STILL  APPEAR  IN  CALLING 
SEQUENCE) 

REMARKS 

ORIGINAL  MATRIX  A  MUST  BE  REAL  SYMMETRIC  (STORAGE  M0DE=1> 
MATRIX  A  CANNOT  BE  IN  THE  SAME  LOCATION  AS  MATRIX  R 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
NONE 
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METHOD 

DIAGONALIZATIQN  METHOD  ORIGINATED  BY  JACOBI  AND  ADAPTED 
BY  VON  NEUMANN  FOR  LARGE  COMPUTERS  AS  FOUND  IN  ' MATHEMAT I CA 
METHODS  FOR  DIGITAL  COMPUTERS'*  EDITED  BY  A.  RALSTON  AND 
H .  S .  WILF  *  JOHN  WILEY  AND  SONS*  NEW  YORK  *  1962*  CHAPTER  7 


SUBROUTINE  EIGEN ( A » R * N * MV > 
REALM  A ( 15  >  *  R ( 15 ) 


IF  A  DOUBLE  PRECISION  VERSION  OF  THIS  ROUTINE  IS  DESIRED  *  THE 
C  IN  COLUMN  1  SHOULD  BE  REMOVED  FROM  THE  DOUBLE  PRECISION 
STATEMENT  WHICH  FOLLOWS. 

DOUBLE  PRECISION  A » R » ANORM . ANRMX . THR » X » Y , SINX * SI NX2 * COSX * 

1  C0SX2 * SINCS *  RANGE »DSQRT*DABS 

THE  C  MUST  ALSO  BE  REMOVED  FROM  DOUBLE  PRECISION  STATEMENTS 
APPEARING  IN  OTHER  ROUTINES  USED  IN  CONJUNCTION  WITH  THIS 
ROUTINE. 

THE  DOUBLE  PRECISION  VERSION  OF  THIS  SUBROUTINE  MUST  ALSO 
CONTAIN  DOUBLE  PRECISION  FORTRAN  FUNCTIONS.  SORT  IN  STATEMENTS 
40*  68*  75*  AND  78  MUST  BE  CHANGED  TO  DSQRT .  ABS  IN  STATEMENT 
62  MUST  BE  CHANGED  TO  DABS.  THE  CONSTANT  IN  STATEMENT  5  SHOULD 
BE  CHANGED  TO  1.0D-12. 


GENERATE  IDENTITY  MATRIX 

5  RANGE*1 . OE-6 

IF(MV-l)  10*25*10 
10  IQ*-N 

DO  20  J=1 *  N 
IQ^IQ+N 
DO  20  1  =  1  *  N 
I J=IQ+I 
R ( I J ) =0 . 0 
IF < I-J)  20*15*20 
15  R(IJ)*1.0 
20  CONTINUE 

COMPUTE  INITIAL  AND  FINAL  NORMS  (ANORM  AND  ANORMX) 

25  ANQRM*0.0 
DO  35  1*1  *  N 
DO  35  J«I*N 
IF(I-J)  30*35*30 
30  IA»I+( J*J-J>/2 

ANQRM“ANORM+A< IA)tA(IA) 

35  CONTINUE 

IF ( ANORM )  165*165*40 
40  ANQRM»1 . 4144SQRT (ANORM) 

ANRMX«ANORM*RANGE/FLOAT(N> 
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INITIALIZE  INDICATORS  AND  COMPUTE  THRESHOLD.  THR 

IND*0 

THR=ANQRM 

45  THR*THR/FLQAT (N) 

50  L  =  1 
55  M=L+ 1 

COMPUTE  SIN  AND  COS 

60  MQ*(M*M-H)/2 
LQ»<L*L-L)/2 
LM*L+MQ 

62  IF (  ABS < A ( LM ) ) -THR )  130.65.65 
65  IND=1 
LL=L+LQ 
MM*M+MQ 

X*0.5*(A(LL)-A(MM)  ) 

68  Y=-A<LM)/  SQRT<A(LM)*A(LM)+X*X> 

IF ( X )  70.75.75 
70  Y=-Y 

75  5INX=Y /  SORT ( 2 ♦ 0*  C 1 . 0+ (  SORT < 1 . 0-Y*Y > > )  ) 
SINX2=SINX*SINX 
78  COSX*  SORT (1 .0-S1NX2) 

C0SX2=C0SX*C0SX 
SINCS  *SINX*CQSX 

ROTATE  L  AND  M  COLUMNS 

ILQ»N*(L-1) 

IMQ=N*(M-1) 

DO  125  I=»1»N 
IQ*(I*I-I)/2 
IF(I-L)  80.115.80 
80  IF(I-M)  85.115.90 
85  IM»I+MQ 
GO  TO  95 
90  IM*M+IQ 

95  IF(I-L)  100.105.105 
100  IL*I+LQ 
GO  TO  110 
105  IL*L+IQ 

110  X*A(IL)»COSX-A(IM)*SINX 

A(IM)»A(lL)«SINX+A(IM)*COSX 
A ( IL ) =X 

115  IF(MV-l)  120.125.120 
120  ILR=ILQ+I 
IMR-IMQ+I 

X=R(ILR)*COSX-R<IMR)#SINX 
R( IMR)*R( ILR)#SINX+R< IMR)*COSX 
R ( ILR ) aX 
125  CONTINUE 

X-2.0*A(LM)*SINCS 

Y=A(LL)*C0SX2+A<MM)*SINX2-X 

X=A<LL>*SINX2+A<MM)*C0SX2+X 

A(LM)=<A<LL)-A<MM> ) *S INCSt A < LM > * ( C0SX2-SINX2 ) 
A<LL>*Y 
A ( MM ) ®X 
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TESTS  FOR  COMPLETION 


TEST  FOR  M  =  LAST  COLUMN 

130  IF(M-N)  135»140tl35 
135  M*M+1 

60  TO  60 

TEST  FOR  L  *  SECOND  FROM  LAST  COLUMN 

140  IF (L-<N-1 ) )  145 » 150» 145 
145  L*L+1 

GO  TO  55 

150  IF(IND-l)  160 » 155 » 160 
155  IND=0 

GO  TO  50 

-  -  COMPARE  THRESHOLD-  I4ITH.  FINAL  -N.QRH  - 
160  IF (THR-ANRMX)  165»165.45 


SORT  EIGENVALUES  AND  EIGENVECTORS 

165  IQ=-N 

DO  185  1*1  »N 
IQ*IQ+N 

LL*I+(I*I-I)/2 

JQ=N*<I-2> 

DO  185  J*I » N 
JQ* JQ+N 

MM=J+< J*J-J)/2 
IF  <  A ( LL ) -A  <  MM ) )  170»1S5»185 
170  X*A<LL) 

A ( LL ) *A  <  MM ) 

A(MM)*X 

IF(MV-l)  175  r 185  > 175 
175  DO  180  K*l t N 
ILR“IQ+K 
IMR=JG+K 
X=R  < ILR ) 

R ( ILR ) *R ( IMR ) 

180  R ( IMR ) *X 
185  CONTINUE 
END 
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BOW  TO  DESIGN  A  NEW  POTENTIAL 


- rh«  potential  Function  used  in“thi«  hfoaram  IS  stored - 

in  the  file  PUT. FOR.  This  particular  Function  arllies  only 
To  system  with  cuo i c  s rmmetrr .  it  must  ITB  modified  for 
performina  calculations  on  other  symmetry  sroups>  but  it 
will  worK  For  the  ootaneoral  case  as  well  as  the  8  of  T2 
Fold  cubic  cases.  The  reason  for  this  is  that  the  lowest 
order  non  spherical  term  Tin  the  potential  due  to  an  octahedral 
coordination  of  ions  is  the  same  as  the  cubic  crystal  field. 

To  determine  what  tne  appropriate  potential  Func t i flh  for 
a  particular  problem  is  consider  the  followina  example. 

Suppose  we  have  an  ion  at  the  poTition  x»-a»  y»Oi  z*o,. This  Ton 
will  contribute  to  the  potential  a  term 


V(x»  y,z)-eC  <a+x)**2+y**2+z**23**-l/2 


which  can  be  expanded  in  a  Taylor  series  to  aive 


V(x, yrz)»e/a( l-x/a-r**2/2a**2+3/2*x**2/a**2+. . . > 


IF  we  place  additional  charaes  at  the  positions  required 
To  Form  an  octahedral  coordination  we  Find  that  a  lot  oF 
cancellation  will  occur,  ‘ihe  First  term  with  nonspherical  symmetry 


x. y » z ) =y»e/ Maees ) * ( x**4*y**e+z**4-3/3*r**4  J 
By  terms  with  spherical  symmetry  we  mean  terms  line 
V<  x » y.z)*x**2+y**2+z**Z 
"or  powers  oF  this. 

This  sort  of  Taylor  series  expansion  of  the  potential 
is  necessary  in  crystal  Field  calculations.  Usina  tnei  Full 
potential  will  not  aive  the  relevant  splittina. 

In  order  to  study  a  particular  distortion  or  other 
symmetry  aroup  the  first  step  is  to  determine  the  symmetry 
aroup.  A  useful  reference  is  TinK ham's  booK  "Group  Theory 

and  Quantum  Mechanics"  . _ 

Once  you  have  the  potential  function  m  it's  laylor  series  form 
the  next  step  is  to  write  a  short  proaram  in  Fortran  to  use 
with  the  other  proarae.  The  proaram  PUTTFUR  on  this  diskette 
should  serve  as  a  model,  lo  aet  a  copy  of  it  type  in 
lYPk  PUl.l-UR. 

Ihe  next  step  is  to  compile*  this  proaram.  The  instructions 
Tn  the  Fortran  manual  should  be  consulted  TF  the  ronowina  directions 

are  not  helpful. _ _ 

T)  Type  in  FORMAT  51 

2)  insert  a  new  disKette  into  the  riaht  hand  disK  drive 

&)*  type  in" COPY"  AI CRYSTAL. OBJ  - 

this  will  copy  the  object  code  for  the  main 
proaram  onto  the  new  disKette 
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4)  t  ype  in  COPY  AILD.UJM  B 

this  will  copy  the  editor 

5)  tyre  in  COPY  A.'POT.FOR  B: 

this  copies  the  potential  Function 

6)  type  in  BlfcU  B : POT .  P  OR 

this  allows  you  to  edit  the  program 
_ refer  to  the  DOS  manual  for  usina  the  editor _ 

7)  once  you  have  edited  your  program  remove  this  diskette  from 

8)  from  the  left  driue  and  insert  the  disk  labelled  FORI _ 

9)  type  in  SI 

10)  type  in  A1F0R1 _  _ _ 

Til  when  the  compiler  ask  you  for  the  name  o?  the  source  code 
_ type  in  POT  _ 

12)  then  hit  the  return  3  additional  times 

13)  when  the  compiler  Finishes  replace  the  disk  on  the  left _ 

with  the  disk  labelled  for2 

14)  type  in  AH-UR2 _ _  _ 

13)  when  this,  finishes  the  second  pass  of  the  compiler 
16)  relace  the  disk  on  the  left  with the  disk  labelled 


a  running  r rosrae  on  the  disk  in  the  right  hand 
drive  called  CRYSTAL. bXh  you  can  verify  this  by 


typing  in  DIR  81 

21)  copy  this  onto  the  disk  on  the  left  by  typing 
COPY  8 1CRYSTAL.EXE  A1 

22)  replace  the  disk  on  the  right  with  the  data  disk 


SOURCE  CODE  FOR  POTENTIAL  FUNCTION 


FUNCTION  POT <  X »  Y  #  Z ) 

COMMON  /PARAM/  ZE,ZO,BD 

r2»k*h+y*y+2*z - 

R4»R2*R2 

T=x*X*)<*X+V*V*Y*Y+2*Z*Z*2-3 . 0*R4/5 . 0 
POl =35.0*ZE/4.0/(BD**5.0)#T 

“ENC 


RANDOM  NUMBER  FILE  GENERATION 


.0  RANDOMIZE _ 

20  OPEN  "RANDOM"  FOR  OUTPUT  AS  #1 

30  FUR  I»1  TO  S003 _ 

40  X*RND 

SO  Y»RND _ _ 

bO  Z»RND 

70  PR1NT#1 t USING" . ####### " ?  X  ?  Y  ?  Z 

bo  prini  using". #######m;x;y;z 

SO  NEXT  1 _ 

100  BEEP :sttp: BEEP: BEEP 

105  CLOSE#  1 _ 

no  end 
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SAMPLE  ROM  OP  PROGRAM  SHOWING  EFFECT  OF  MONTE  CARLO 
POINTS  ON  CONVERGENCE  OF  EIGENVALUES 

cat  t ink 1000 

ENTER  04  FOR  TETRAHEDRAL  COORDINATION 
ENTER  06  FOR  OCTAHEDRAL  COORDINATION 
ENTER  08  FOR  8-FOLD  CUBIC  COORDINATION 
ENTER  12  FOR  12-FOLD  CUBIC  COORDINATION 
ENTER  0  TO  RUN  WITH  INTERNAL  PARAMETERS 
ENTER  1  TO  SPECIFY  NEW  PARAMETERS 
POSITIONING  IONS  IN  THE  LATTICE 
INITIALIZING  CHARGE  ON  OXYGEN  IONS 
HOW  MANY  MONTE-CARLO  POINTS  (15)  LESS  THAN  32000 
COMPUTING  MATRIX  ELEMENT 
-0 . 75801 Id  00 
COMPUTING  MATRIX  ELEMENT 
-0. 193706d-01 
COMPUTING  MATRIX  ELEMENT 
-0.675547d  00 
COMPUTING  MATRIX  ELEMENT 
0 . 161519d-01 

COMPUTING  MATRIX  ELEMENT 
0 . 328048d-01 

COMPUTING  MATRIX  ELEMENT 
-0 . 581051d  00 
COMPUTING  MATRIX  ELEMENT 
-0. 189928d-01 
COMPUTING  MATRIX  ELEMENT 
0.284910d-01 

COMPUTING  MATRIX  ELEMENT 
0.374327d-02 

COMPUTING  MATRIX  ELEMENT 
0 . 828496d  00 

COMPUTING  MATRIX  ELEMENT 
-0.782377d-01 
COMPUTING  MATRIX  ELEMENT 
-0  « 683282d-01 
COMPUTING  MATRIX  ELEMENT 
-0.457849d-01 
COMPUTING  MATRIX  ELEMENT 
0  .  743769d-01 

COMPUTING  MATRIX  ELEMENT 
0.903533d  00 
COMPUTING  EIGENVALUES 
THE  CALCULATION  PARAMETERS  ARE . 

CHARGE  ON  OXYGEN  2.00 
EFFECTIVE  NUCLEAR  CHARGE  2.00 
LATTICE  PARAMETER  2.00 
THE  EIGENVALUES  ARE . 

0 . 954882d  00  0.786077d  00  -0.573409d  00  -0.680228d  00  -0.769904d  00 

6 

1000 

.00 

.00 

.00 

0 . 954882d  00  0.786077d  00  -0.573409d  00  -0.680228d  00  ~0.769904d  00 
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cat  tink4000 

ENTER  04  FOR  TETRAHEDRAL  COORDINATION 
ENTER  06  FOR  OCTAHEDRAL  COORDINATION 
ENTER  08  FOR  8-FOLD  CUBIC  COORDINATION 
ENTER  12  FOR  12-FOLD  CUBIC  COORDINATION 
ENTER  0  TO  RUN  WITH  INTERNAL  PARAMETERS 
ENTER  1  TO  SPECIFY  NEW  PARAHETERS 
POSITIONING  IONS  IN  THE  LATTICE 
INITIALIZING  CHARGE  ON  OXYGEN  IONS 
HOW  MANY  MONTE-CARLO  POINTS  <I5)  LESS  THAN  32000 
COMPUTING  MATRIX  ELEMENT 
-0.684433d  00 
COMPUTING  MATRIX  ELEMENT 
0 . 622392d-02 

COMPUTING  MATRIX  ELEMENT 
-0.653617d  00 
COMPUTING  MATRIX  ELEMENT 
0. 115858d-01 

COMPUTING  MATRIX  ELEMENT 
0. 195847d-01 

COMPUTING  MATRIX  ELEMENT 
-0 . 624519d  00 
COMPUTING  MATRIX  ELEMENT 
-0. 130734d-01 
COMPUTING  MATRIX  ELEMENT 
0 . 993926d-02 

COMPUTING  MATRIX  ELEMENT 
0 . 152267d-01 

COMPUTING  MATRIX  ELEMENT 
0 . 81 4384d  00 

COMPUTING  MATRIX  ELEMENT 
-0 . 482886d-01 
COMPUTING  MATRIX  ELEMENT 
0 . 171561d-01 

COMPUTING  MATRIX  ELEMENT 
0 • 381 198d-01 

COMPUTING  MATRIX  ELEMENT 
0 . 217708d-01 

COMPUTING  MATRIX  ELEMENT 
0 . 748680d  00 
COMPUTING  EIGENVALUES 
THE  CALCULATION  PARAMETERS  ARE . 

CHARGE  ON  OXYGEN  2.00 
EFFECTIVE  NUCLEAR  CHARGE  2.00 
LATTICE  PARAMETER  2.00 
THE  EIGENVALUES  ARE . 

0 . 822025d  00  0.744228d  00  -0.613175d  00  -0,663467d  00  -0.689116d  00 

6 

4000 

.00 

.00 

.00 

0 . 822025d  00  0.744228d  00  -0.613175d  00  -0.663467d  00  -0.689116d  00 
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cat  tink.32000 

ENTER  04  FOR  TETRAHEDRAL  COORDINATION 
ENTER  06  FOR  OCTAHEDRAL  COORDINATION 
ENTER  08  FOR  8-FOLD  CUBIC  COORDINATION 
ENTER  12  FOR  12-FOLD  CUBIC  COORDINATION 
ENTER  0  TO  RUN  WITH  INTERNAL  PARAMETERS 
ENTER  1  TO  SPECIFY  NEW  PARAMETERS 
POSITIONING  IONS  IN  THE  LATTICE 
INITIALIZING  CHARGE  ON  OXYGEN  IONS 
HOW  MANY  MONTE-CARLO  POINTS  (15)  LESS  THAN  32000 
COMPUTING  MATRIX  ELEMENT 
-0.614619d  00 
COMPUTING  MATRIX  ELEMENT 
-0. 197246d-02 
COMPUTING  MATRIX  ELEMENT 
-0 » 627631d  00 
COMPUTING  MATRIX  ELEMENT 
-0. 187123d-01 
COMPUTING  MATRIX  ELEMENT 
-0. 181665d-01 
COMPUTING  MATRIX  ELEMENT 
-0 . 61 1 702d  00 
COMPUTING  MATRIX  ELEMENT 
-0 . 124754d-02 
COMPUTING  MATRIX  ELEMENT 
0.927798d-02 

COMPUTING  MATRIX  ELEMENT 
-0 . 95297 Id- 02 
COMPUTING  MATRIX  ELEMENT 
0.801265d  00 

COMPUTING  MATRIX  ELEMENT 
0 • 578322d-03 

COMPUTING  MATRIX  ELEMENT 
0 . 122425d-0 1 

COMPUTING  MATRIX  ELEMENT 
-0.224950d-02 
COMPUTING  MATRIX  ELEMENT 
-0 • 619347d-02 
COMPUTING  MATRIX  ELEMENT 
0.788407d  00 
COMPUTING  EIGENVALUES 
THE  CALCULATION  PARAMETERS  ARE . 

CHARGE  ON  OXYGEN  2.00 
EFFECTIVE  NUCLEAR  CHARGE  2.00 
LATTICE  PARAMETER  2.00 
THE  EIGENVALUES  ARE.,.,. 

0 . 803822d  00  0.786089d  00  -0.590134d  00  -0.619015d  00  -0.645042d  00 

6 

2000 

.00 

.00 

.00 

0 . 803822d  00  0.786089d  00  -0.590134d  00  -0.619015d  00  -0.645042d  00 
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